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The microscopic magnetic properties of high-quality GdN thin films have been investigated by 
electron spin resonance (ESR) and ferromagnetic resonance (FMR) measurements. Detailed temper¬ 
ature dependence ESR measurements have shown the existence of two ferromagnetic components 
at lower temperatures which was not clear from the previous magnetization measurements. The 
temperature, where the resonance shift occurs for the major ferromagnetic component, seems to 
be consistent with the Curie temperature obtained from the previous magnetization measurement. 

On the other hand, the divergence of line width is observed around 57 K for the minor ferromag¬ 
netic component. The magnetic anisotropies of GdN thin films have been obtained by the analysis of 
FMR angular dependence observed at 4.2 K. Combining the X-ray diffraction results, the correlation 
between the magnetic anisotropies and the lattice constants is discussed. 


2 


I. INTRODUCTION 


Ferromagnetic semiconductors are expected to be a key material for the future spintronics.-^ GdN is one of these 
ferromagnetic semiconductors and it is particularly interesting due to its partially filled 4f and 5d orbitals with 
saturation moment of 7 ^b/ Gd 3+ £ Therefore, GdN has been the object of a series of theoretical and experimental 
studies since more than a half century*^— 

However, it is well known from a number of studies of bulk GdN in the 60’s and 70’s^A— that it is very difficult to 
obtain the high-quality bulk GdN because nitrogen vacancies and oxygen can damage it very easily. For example, it 
was reported that there is a strong decrease in the magnetic moment of bulk GdN even with few percent of oxygen— 
The decrease in Curie temperature of bulk GdN was confirmed with a range of oxygen concentration* 1 -- Cutler et al 
also reported that the nitrogen vacancies make the hysteresis effect and the remenance much smaller.— Understanding 
the properties of ’’pure” GdN is still challenging because of the difficulty to produce high quality single crystals. 

However, the situation has changed owing to the advanced technology of the thin film synthesis since the 2000s ; 7 i 13 i 22 
These thin films are capped by AIN on GdN to make surface smooth and to restrict the oxygen contamination. These 
early studies of GdN thin film reported that the properties of the GdN thin film, contrary to bulk one, are very 
sensitive to their epitaxial strain, structural distortion and surface effect for nanocrystalline films. From the view 
point of the development in spintronics devices, therefore, it is very important to understand these parameters and 
how they affect the magnetic and optical properties of the high-quality GdN thin film. 

More recently, H. Yoshitomi et al. and R. Vidyasagar et al. have studied the optical and magnetic properties 
in epitaxial AIN/GdN/AIN double heterostructures grown by reactive radio-frequency (rf) sputtering under ultra- 
pure conditions— r— For example, their high-quality GdN thin film of 95 nm showed the indirect and direct optical 
transitions, and the considerable size effects of the optical band gap were observed with a decrease in the GdN 
thickness. They also investigated the saturation magnetization and Curie temperature estimated by Arrott plots as a 
function of the thickness of GdN ; 13 ' 23 However, few cases, except for the ferromagnetic resonance (FMR) measurement 
by K. Khazen et al?%, studied the microscopic magnetic properties of these GdN thin films. 

In this study, we investigate the microscopic magnetic properties of high-quality GdN thin films by the detailed 
temperature dependence of electron spin resonance (ESR), and the angular dependence of FMR at 4.2 K. 


II. EXPERIMENTAL DETAILS 


We investigate the micro magnetic properties of three GdN samples. One has the thickness of 95nm whose optical 
and macro magnetic properties has been investigated by H. Yoshitomi et al22- We call this sample “08GdN” in this 
paper. The other samples has the thicknesses of 29 nm and 97 nm, respectively; we call these two samples “lOGdN” 
in this paper. All samples were grown on c-sapphire (0001) substrates at 500 °C by reactive radio-frequency (RF) 
magnetron sputtering^ in an ultrahigh vacuum chamber. The input RF power was 250 W. AIN/GdN/AIN double 
heterostructures were used to avoid oxidation— The growth chamber equipped with multitargets for AIN and GdN was 
separated from the substrate introduction chamber to avoid oxidation of the target when introducing the substrate. 
Al(99.99%) and Gd(99.9%) were used as metal targets. We used an ultrapure (99.9999%) gas mixture of argon and 
nitrogen for reactive growth. For the synthesis of 08GdN thin film, the partial pressure ratio of argon and nitrogen 
was even, and the total sputtering pressure was 5 Pa. For the syntheses of lOGdN thin films, on the other hand, the 
partial pressure ration was 9:6, and the total sputtering pressure was 6 Pa for the purpose to decrease the number 
of nitrogen vacancies. However, the transmission and absorption spectral measurement showed that the number of 
free carrier in lOGdN is more than that in 08GdN. It’s not known exactly why the number of nitrogen vacancy in 
lOGdN is more than that in 08GdN— The X-ray diffraction measurement showed that the lattice constant along 
(111) direction for 08GdN sample is a = 0.507 nm. The 29nm and 97nm thin films of lOGdN have a = 0.506 nm 
and a = 0.507 nm, respectively. These lattice constants are longer than the bulk value a = 0.4998 nm. In addition, 
we also confirmed that the lattice constant along (200) direction is smaller than the bulk value a = 0.250 nm; for 
example, the value for 95 nm thickness of 08GdN is a = 0.249 nm. Therefore, our GdN thin films have uniaxial lattice 
distortion. 

Our ESR/FMR measurements were performed by the Bruker X-band ESR spectrometer EXM081 at Center for 
Supports to Research and Education Activities, Kobe University, with 100kHz field modulation using a TE103 rect¬ 
angular cavity in the temperature range of 4.2 K to 300 K. We show the geometry of the FMR measurements in 
Fig. 1, which is the same condition as that in Khazen’s paper* 22 The GdN samples lie in the x-z plane, and the y-axis 
is parallel to the growth face direction [111] of our thin films. The out-of plane variation of the external magnetic 
field is in the xy plane. 
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Figure 1. (Color online) The geometry of the FMR measurements. M and H are spontaneous magnetization and external 
magnetic field, respectively. 6(6h), ifi(ifiH) are polar and azimuthal angles for M(H) vector. 


III. EXPERIMENTAL RESULTS 

Firstly, we show the temperature dependence of the resonance field. We successfully obtained the result of Fig. 2 
(a) from the spectrum fitting, using the following double lorentzian equation with respect to the ESR spectrum at 
each temperature. Here, we would like to point out that the observed ESR signal is the differential curve. Total 
intensity I as a function of external magnetic field H is written by 

T(m = /^[(g-gQ/CAgpp!^)] 

1 1 [3 + {(H — Hi) / (AHppi/2)} 2 ] 2 

I' m2 [(H-H 2 )/{AH P p 2 /2)] 

[3 + {(H — H 2 ) / (AHpp2/2)} 2 ] 2 

+m + m'H, (1) 

where I' mi is the intensity, Hi is the resonance field, AHppi is the line width for component i and m and in' are 
background. These parameters are determined by the fitting experimental data. As an example, we show the ESR 
spectrum and the fitting result at T =40 K for 95 nm thickness of 08GdN in the inset of Fig. 2 (a). Note that 
the external magnetic held was applied to in the plane (ipn = 0). Fig. 2 (a) shows the temperature dependence 
of resonance held for 95 nm thickness of 08GdN and 97nm thickness of lOGdN. We estimate the <?-factor g ~ 1.96 
by using the resonance held value 3524.9 G at the highest temperature 260 K. This g value is consistent to that of 
Lande g-factor gL = 2 of Gd 3+ whose total orbital angular momentum L is 0. At the low temperature region, we 
can see clearly two kinds of phases not only for 08GdN but also for lOGdN which has the larger number of nitrogen 
vacancy. Therefore, the origin of the phase separation is not coming from the nitrogen vacancy. We have already 
confirmed the existence of such two kinds of phases in the other GdN samples^ The Curie temperature (T c ) for the 
95 nm thickness of 08GdN has been reported about 37 K by using Arrott plot analysis^, therefore, the shift in the 
resonance field around 40 K (res. 2 in Fig. 2 (a)) comes from a dominant part of the magnetization of GdN thin 
film in the ferromagnetic phase. In the present study, we also confirmed that the T c value for the 97 nm thicknesses 
of lOGdN is about 29 K by using Arrott plot analysis. The difference of the T c between 08GdN and lOGdN comes 
from the number of nitrogen vacancy^ Careful observation of Fig. 2 (a) tells us that the resonance shift of 10 GdN 
(res. 2) begins at lower temperature than that of 08 GdN, and which is consistent that the Curie temperature for 10 
GdN is lower than that of 08 GdN. Therefore, the shift in the resonance held around 30 K (res. 2) for lOGdN also 
comes from a dominant part of the magnetization in the ferromagnetic phase. The second shift at higher temperature 
side (^70K (res. 1)) is originated from another ferromagnetic phase, which cannot be ascribed to the short-range 
correlation of spins at T > T c because without the phase separation we cannot observe two ESR’s in the intermediate 
temperature region (25^57 K). The high-Tc phase may come from the interface because the contribution of res. 1 to 
the static magnetization is less than 1% for 08GdN at 50 K where the shift of res. 1 is close to the saturation while 
the shift of res. 2 has just started^ However it requires further investigation. Here we would like to emphasize that 
the observation of the two kinds of phases as in the case of Fig. 2 (a) suggests ESR measurements can easily detect 
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Figure 2. (Color online) (a) Temperature dependences of the resonance field for 95 nm thickness of 08GdN and 97 nm thickness 
of lOGdN. The word of “res. 1(2)” means the resonance shift at higher (lower) temperature region. The inset is the fitting 
result at T = 40 K for 95 nm thickness of 08 GdN by using double Lorentzian equation (1). (b) and (c) are temperature 
dependences of the line width and of integrated intensity for 95 nm thickness of 08GdN, respectively. 


the microscopic properties such as the phase separation which is difficult to be observed by macro measurements. 
We also note that it was difficult to measure the resonance field for 10 GdN samples at higher temperature region 
because the nitrogen vacancy provides carriers at higher temperature which causes the decrease of Q-factor in ESR 
measurements. 

Next, we investigated the temperature dependence of the line width for 95 nm thickness of 08GdN. Fig. 2(b) shows 
the results about the line width. Owing to ESR measurements, we confirm the decreasing behaviors of the line 
width AHppi and AHpp 2 below 57 K which corresponds to the mid-point of res. 1 resonance shift in Fig. 2(a). The 
decreasing behavior with decreasing temperature in the ferromagnetic phase is well known as a typical property in 
the ferromagnetic region^ It is also interesting that the divergence behavior of AHppi can be observed clearly near 
at 57 K. This divergence of the line width indicates the presence of spin fluctuations near T c .— The similar increase 
behavior near above T c has been also confirmed roughly in typical FM thin films.* 31 However, no such divergence 
behavior is observed for AHpp 2 suggesting the different spin dynamics for res. 1 and res. 2. We also investigated 
the temperature dependence of the integrated intensity as shown in Fig. 2(c). Here, integrated intensity Si for each 
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component i can be calculated by the line width Affpp,; and intensity I' mi . It is well known in ESR that 5) can be 
estimated by I' mi (AHppi ) 2 . According to this result, we can determine which of two separated phases is smoothly 
connected to the dominant part in the ferromagnetic ground state. Therefore, we can say that “res. 2” in Fig. 2 (a) 
and AHpp 2 in Fig. 2(b) correspond to the dominant part in ferromagnetic state. 

Next, we investigated the angular dependence of FMR. For example, Fig. 3 shows the angular dependence of the 
FMR for the 08GdN at 4.2K where is varied from —90° to 90°. Here, we measured the FMR in the interval of 
10 degrees. We confirm that the resonance field is very sensitive to the angular variation. Applying the equation (1) 
to analyze data in Fig. 3 we get the angular dependence of the line width in Fig. 4. Owing to the closely-spaced 
angulars we took, we see the peak structures clearly at ipn ~ 90° in Fig. 4. We also gain the angular dependence of 
the resonance field. Fig. 5 presents the angular dependence of the resonance field for our three samples obtained by 
the Lorentzian fitting. In order to investigate the film thickness dependence, we add the data for 29 nm thickness of 
lOGdN sample. All sample shows that the resonance field was maximized when external magnetic field was applied 
in the direction of out-of-plane, and minimized when in the direction of in-plane. This result consistent with typical 
FMR spectral, and with behaviors of Khazen’s samples^. A careful observation of Fig. 5 enable us to confirm the 
resonance field value for 29 nm thickness of lOGdN is lower than those for the other samples. This behavior means 
that the magnetic anisotropy for 29 nm thickness of lOGdN is different from the other samples. 

Finally, we analyzed the magnetic anisotropy for our GdN thin film from the angular dependence of the FMR 
spectra. We use “Smit-Bcljer formalism”— for our analysis which is applicable to thin film of cubic symmetry 
allowing for a possible uniaxial deformation, because our GdN samples also have uniaxial anisotropy in the process 
of synthesis. In this case, the energy density E can be written by 


E = -MH + K^alal + a\a\ + a 2 3 aj) + (2t tM 2 - K n )a 2 2 , (2) 

which represent the Zeeman interaction, the magnetic anisotropic energy, and the demagnetization energy. Here, 
Ki is the fourth order cubic magnetocrystalline anisotropy constant, and K u is the second order uniaxial anisotropy 
constant. o.i is the direction cosines of the magnetization M relative to the cubic crystal axes, and H is the applied 
field (See Fig. 1). In order to analyze the results in Fig. 5, we can generally fix 6 = Oh = 7r/2 in Fig. 1 and use the 
following equations. One is the static equilibrium orientation of the magnetization 

Hsin(<p H -ip) = (4? tM - ^)sin(^)cos (ip) + |^sin(4<£), ( 3 ) 

and the other is the resonance field equation 


These two fitting equations for 
Beljers equatons^): 


K 2 K 

= [Hcos{ip H -ip) + -^(2 - sin 2 2<^) - (4? tM - -^)sin 2 <p] 

2 o ft 

x [Hcos(iph — tp) + -^-cos4y> + ( 4?r -^--^-)cos2(^]. (4) 

FMR measurements are derived from the following three resonance conditions (Smit- 
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In Fig. 5 the example of the fitting result for 08GdN sample is shown by the black line. The fitting is rather successfull. 
Small deviations between the data and the fitting close to ±90° may be due to the subtle misalignment ot the sample 
to the applied magnetic field. We also performed the same fitting to the obtained data for lOGdN 29 nm and 97 nm 
samples where the fitting lines are not shown in Fig. 5 to avoid the complication in Fig. 5. However, the fittings are 
also rather successful. Table 1 shows our analysis results for magnetic anisotropy constants. Here the magnetization 
M for each sample is obtained from the paper* 2 ' These crystal anisotropies K\ and K u are much different from 
the Khazen’s results^. Our K\ value is almost one third, and K u is two or three times of each value of Khazen’s 
bulk sample, respectively. The reason of the difference comes from the difference in the crystal growth process. The 
Khazen’s GdN samples were deposited on (100) oriented Si substrate and these films were polycrystalline. On the 
other hand, our samples were grown along to the (111) direction of GdN on c-sapphire (0001) substrates by reactive 
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Figure 3. Angular dependence of FMR measurement at T =4.2 K for 95 nrn thickness of 08GdN. 


radio-frequency magnetron sputtering in an ultrahigh vacuum chamber. More concretely, the lattice constant of the 
Khazen’s extended film was increased 2.4% uniformly not along to a specific direction, and the K\ value of 2.4 % 
increased samples are larger than that of Khazen’s bulk sample. The lattice constant along (200) direction of our 
samples were decreased although the lattice constant along (111) direction were increased. In fact, the lattice constant 
along (111) direction for the samples of the 95 nm thickness of 08GdN and 97 nm thickness of lOGdN is is a =0.507 
nm which is larger than the reported bulk value a=0.4998 nm.The lattice constant along (200) direction for 97 nm 
thickness of 08GdN is a=0.249 nm and this value is smaller than the bulk value of a =0.250 nm. In addition, the 
coefficient of thermal expansion of our substrate AIN is 4.4 xl0~ 6 /\ _1 which is larger than that of the Si substrates, 
2.4 xl0 -6 /\ -1 . Therefore, we can naturally accept the difference between our and Khanzen’s samples. We also note 
that the K\ value for 29 nm thickness of lOGdN is larger than those of our other samples. The largeness comes from 
the small lattice constant for 29 nm thickness of lOGdN along (111) direction, a =0.506 nm. In other words, the 
K\ values has the tendency to come close to the bulk value when the lattice constant approaches to the bulk value 
a =0.4998 nm. We should be careful to the fact that these lattice constant values measured by X-ray diffraction are 
just average values. Therefore, we can not discuss about the second order uniaxial anisotropy K u from the view point 
of the lattice constant because the K u values are mainly affected by the interfacial surface of crystal. It is naturally 
expected that our K u values are more sensitive than Khazen’s sample because Khazen’s films are polycrystalline. This 
is the origin that our K u values are larger than those of Khazen’s. The K u value of 95 nm thickness for 08 GdN is 
slightly larger than that of 97 nm thickness for lOGdN. We speculate that it may come from the nitrogen vacancy, that 
means, the strain at the interfacial surface for lOGdN samples was relaxed by the large number of nitrogen vacancy. 
We also comment about the characteristic which our thinner sample of 29 nm thickness for 10 GdN has the largest 
K u value. It is characteristic of ferromagnetic thin films that the thinner thickness sample has the larger value of K u . 
This behavior is well known theoretically and experimentally, for example, Fe/MgO multilayered films^ According 
to these obtained results, the FMR analysis is very useful to obtain the microscopic properties. 
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Figure 4. Angular dependence of line width at T =4.2 K for 95 nm thickness of 08GdN. 



Figure 5. (Color online) Angular dependences of resonance field of FMR for our all samples at T = 4.2 K. The external 
magnetic field was applied out-of-plane at —90° and 90°. The black line is the fitting result for 08GdN sample. 



47 tM (Oe) 

2 K u /M (Oe) 

2 Kx/M (Oe) 

K u (erg/cm 3 ) 

A'l (erg/cm 3 ) 

T(K) 


08GdN 95nm 

24167 

11167 

404 

1.07xl0 7 

3.88x10 s 

4.2 

in this paper 

lOGdN 29nm 

28660 

14760 

394 

1.68xl0 7 

4.49x10 s 

4.2 

in this paper 

lOGdN 97nm 

24027 

10527 

380 

l.OlxlO 7 

3.63x10 s 

4.2 

in this paper 

Bulk film 

22220 

5759 

1292 

5.09x10® 

1.14x10® 

4.0 

ref. 23 

2.4 % extended film 

15620 

2897 

2252 

1.8 xlO® 

1.4 x 10® 

4.0 

ref. 23 


Table I. Spontaneous magnetization M, the fourth order cubic magnetocrystalline anisotropy constant K\ and the second order 
uniaxial anisotropy constant K u for our samples and Khazen’s. 


IV. CONCLUSION 


We have investigated microscopic properties of high-quality GdN thin films. Detailed temperature dependence 
ESR measurements have been performed for the first time and they showed the existence of two ferromagnetic 
components at lower temperatures. It also showed that the temperature, where the resonance shift occurs for the 
major ferromagnetic component, seems to be consistent with the Curie temperature obtained from the previous 
magnetization measurement. On the other hand, the divergence of line width is observed around 57 K for the minor 
ferromagnetic component. We have also determined the fourth order cubic magnetocrystalline and second order 
uniaxial anisotropies of our GdN samples from the angular dependence of FMR measurements observed at 4.2 K. Our 
analysis by Smit-Beljer formalism have clarified that the cubic anisotropy is very sensitive to the lattice constant of 

















thin film and the uniaxial anisotropy values depend on the thickness of thin films strongly. 
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